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Abstract 

The theory of weakly bound cluster breakup, like halo nucleus 
breakup, needs an accurate treatment of the transitions from bound 
to continuum states induced by the nuclear and Coulomb potentials. 
When the transition probability is not very small, a non-perturbative 
framework might be necessary. Nuclear excitation dominates at small 
impact parameters whereas the Coulomb potential being long range 
acts over a larger impact parameter interval. In this article, we propose 
an effective breakup amplitude which meets a number of requirements 
necessary for an accurate quantitative description of the breakup re- 
action mechanism. Furthermore our treatment gives some insight on 
the interplay between time dependent perturbation theory and sud- 
den approximation and it allows to include the nuclear and Coulomb 
potentials to all orders within an eikonal-like framework. 

1 Introduction 

Break up of halo nuclei is a stimulating field for reaction mechanism studies 
where accurate reaction theories are needed in order to extract spectroscopic 



1 



information from experimental data. These theories should treat together 
very different excitations induced by electromagnetic and nuclear fields to 
all orders. Recently a number of interesting papers have appeared in which 
the problem of nuclear and Coulomb breakup is solved numerically either 
by a direct solution of the single particle Schrodinger equation [T]- [7], or by 
DWBA type of approaches [8J or by an approximate treatment of the cou- 
pled equations in the continuum P . Of these papers only Refs. [21 El El E] 
have treated at the same time the nuclear and Coulomb processes. Still one 
needs analytical models to complete the understanding of the mechanisms 
involved and test approximations which could help reduce the computational 
difficulties and help extending the models to more structured clusters. Var- 
ious analytical solutions of the Coulomb breakup problem exist where the 
problem of the higher order effects has been studied pQ. In this paper, we 
propose an effective amplitude for breakup to all orders in the interactions 
which reveals the interplay between sudden and time dependent perturbation 
theory. 

Nuclear breakup is a short time process, several observables of which are 
reasonably well described by the sudden approximation QUI . In the case 
of the electromagnetic field, things are more complicated since the Coulomb 
potential is a long range field. Hence, in Section 2 of this paper we present 
the theoretical framework for the simultaneous treatment of nuclear and 
Coulomb interactions. The effective amplitude we propose is introduced in 
Section 3 where we show also how one can establish the accuracy of different 
approximations based on the relative behavior of two parameters. The last 
part of this article is dedicated to the discussion of some calculations and to 
their comparison with experimental data. 

2 Eikonal theory of nuclear and Coulomb breakup 

In a recent paper JT] we presented a full description of the treatment of the 
scattering equation for a projectile which decays by single neutron breakup 
due to its interaction with the target. There it was shown that within the 
semiclassical approach for the projectile-target relative motion, the amplitude 
for a transition from a nucleon bound state ipi in the projectile to a final 
continuum state i/jf is given by 




(1) 
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where V is the interaction responsible for the neutron transition to the con- 
tinuum. 

For light targets the recoil effect due to the projectile-target Coulomb 
potential is rather small and the interaction responsible for the reaction is 
mainly the neutron-target nuclear potential. In the case of heavy targets the 
dominant reaction mechanism is Coulomb breakup. The Coulomb force does 
not act directly on the neutron but it affects it only indirectly by causing the 
recoil of the core. Therefore the neutron is subject to an effective force which 
gives rise to an effective Coulomb dipole potential V e s(r, R(£)) (cf. Eq. (|23~j) 
of the Appendix). R(£) = ~kd + zvt is the core-target relative distance and 
r is the neutron-core coordinate. In ref . jl lj it was shown that the combined 
effect of the nuclear and Coulomb interactions to all orders can be taken into 
account by using the potential V = V n t+V e {[ sum of the neutron-target optical 
potential and the effective Coulomb dipole potential. If for the neutron final 
continuum wave function we take a distorted wave of the eikonal-type, then 
the amplitude in the projectile reference frame becomes : 



where 0, is the time independent part of the neutron initial wave function 
and i = (l,m) stands for the angular momentum quantum numbers, to = 
{s'f — £oJ A and Eq is the neutron initial bound state energy while s'f is the 
final neutron-core continuum energy, k = (k x , k y , k z ) is a real vector and the 
eikonal phase shift is simply 



where v is the relative motion velocity at the distance of closest approach. 
Integrating by parts Eq. (J2J) leads to the equivalent expression for the breakup 
amplitude : 



Eq. (J2J is appropriate to calculate the coincidence cross section A p — > (A p — 
1) +n. Finally the differential probability with respect to the neutron energy 
and angles can be written as 






(3) 
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where is given by Eq. (J2J and we have averaged over the neutron initial 
state. 

Eq. (J2J) can be in principle an useful alternative to full numerical solutions 
of the Schrodinger equation. In fact it contains all partial waves in the 
final eikonal-like wave function and still the full time dependence, while the 
numerical solutions so far available are often restricted to the first few partial 
waves in the development of the final continuum wave function. Finally the 
differential breakup cross section is given by an integration over core-target 
impact parameters 



and C 2 S is the spectroscopic factor for the initial single particle orbital. 
The effects associated with the core-target interaction have been included by 
multiplying the breakup probability by P c t(d) = \S ct \ 2 [12] the probability 
for the core to be left in its ground state, defined in terms of a core-target 
S-matrix function of d, the core-target distance of closest approach. A simple 
parameterization is P c t(d) = exp(— (In 2) exp[(i? s — d)/a\), where the strong 
absorption radius R s ~ 1A(A^ 3 +A] / ' 3 ) fm is defined as the distance of closest 
approach for a trajectory that is 50% absorbed from the elastic channel and 
a = 0.6 fm is a diffusness parameter. 

There have been already in the literature a large number of papers deal- 
ing with the problem of higher order effects in halo breakup and therefore 
it is important to understand the relation between our model and other ap- 
proaches. By using a first order time dependent amplitude in Eq. we are 
assuming that breakup is a one step process in which the neutron is emitted 
in the continuum by a single interaction with the nuclear target potential 
and by core recoil. The nuclear and Coulomb potential are seen as final state 
interactions which distort the simple plane wave which otherwise would be 
the final continuum state of the neutron. Because of the long tail of the halo 
wave function the overlap between ifi(r) and the potential V(r) is large and 
the potential needs to be treated to all orders. This approach is fully consis- 
tent with the usual treatment of higher order effect by the electromagnetic 





field p. 
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3 Approximation scheme : BBM 1 and BBM 
2 amplitudes 



Since it is already well established that the nuclear potential needs to be 
considered to all orders for weakly bound projectiles, in [11. we studied nu- 
merically only the limits of pure nuclear breakup to all orders, of Coulomb 
breakup to first order and of the coupling between Coulomb to first order 
and nuclear to all orders. On the other hand we argued that the question 
of if and when the Coulomb potential needs to be treated to all orders was 
still under investigation. Here we report on new calculations that we have 
recently performed by using Eq. (J2J) to treat in detail Coulomb higher order 
effects. 

Expanding the time dependent perturbation theory amplitude Eq. (J2J) in 
powers of the eikonal phase shift Eq. 

9cn = £cn + 9 2 cn + 9cn + ■■■ ( 6 ) 

we get a series of partial amplitudes. From here on we avoid the indices 
fi to simplify the notation. Treating separately the nuclear and Coulomb 
potential Eq. (0) reduces to 

= 9&* 1 + S%* 2 + <%*' + .. . (7) 
= gg* 1 + 9$*' + !&*' + .. . (8) 

On the other hand if we make the sudden approximation u = then the 
analogous amplitudes are 

sudd sudd 1 i sudd 2 i sudd 3 i 

9c ~ 9c 9q t yc T • • • {i>) 

sudd suddl , sudd 2 , sudd 3 , f~in\ 

9n — 9n +9n -r • • • yi-v) 

In Eq. (JDJ) each term corresponds to the n th term of the standard eikonal 
approximation to the theory of Coulomb excitations while in Eq. (|7jl 
gQ Ttl is the standard first order perturbation theory. Also we know already 
that the sudden approximation gives an accurate framework for the nuclear 
breakup, then 1 — g^ ddl ~ for each order, such that g^ dd ~ g^ ert - 

Eq. (JSJ) is much easier to calculate than Eq. ((7|), however it has the very 
well known drawback that the first order term leads to a logarithmic diver- 
gence when used in the integral over impact parameters Eq. (|5|i. Under the 
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hypothesis that higher order terms are accurately calculated by the sudden 
approximation, we propose the use of an effective amplitude defined as: 

^BBMl^^sudd_^uddl) + ^e r tl (n) 

We want to treat the nuclear process at the same level of approximation, 
hence, it is straightforward to show that the generalization of g^ BM1 which 
includes also the nuclear potential to all orders and the coupling between the 
nuclear and the Coulomb effective potential is simply 



BBM 1 

ycN 




In the above equation the choice to treat terms higher than the first order 
within the sudden approximation might seem somewhat arbitrary. Then, we 
define another amplitude, BBM 2, for which we keep the time dependence in 
both first and second order terms, sudden approximation being used starting 
from the third order term : 

BBM 2 _ pert 1 . pert 2 . / sudd sudd 1 sudd2\ 

ycN = 9c +3c t yfJcN ~ yc ~ yc ) 

+ (yc POrtl -yc ddl )^ uddl (13) 

This approximation scheme solves several problems encountered in the 
treatment of halo breakup: the already mentioned logarithmic divergence in 
the impact parameter integral due to the the first order sudden approxima- 
tion, the requirement to treat the Coulomb and the nuclear field at the same 
level of approximation and to all orders for small impact parameters and the 
need to use time dependent perturbation theory for large impact parameters 
where the Coulomb field is effective for a long time. A quantitative justifica- 
tion of our approximation scheme is given by the discussion of Fig. ([TJ in the 
next section. On the other hand the discussion of Fig. (2) where the results 
obtained with BBM1 and BBM2 are compared will clarify the treatment of 
higher order terms by the eikonal approximation. 

4 Time dependent framework and its sudden 
approximation 

In this section we give some explicit expressions for the amplitudes discussed 
in Sec. (3). We start by considering the Coulomb term only and in particular 
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the first order approximation for it, thus exp(— iXes) = 1 and V nt = 0, but 
the ujt term is kept in Eq. (@J) (this is the standard first order time dependent 
perturbation theory amplitude) 

#g ertl (k, d;lm) = Q {^wKi (zu) + ivuK (zu) -^j <fi im (k) . (14) 

Here Q = ZfiiZpZxe 2 /hvd is the classical Coulomb momentum transfer to 
the neutron due to the core recoil. Ki and Kq are the usual modified Bessel 

functions. The adiabaticity parameter zu = £f n £ ° d represents the ratio of 
the collision time {d/v) over the nuclear interaction time. If the reaction 
mechanism is such that w is small, then the nuclear interaction time is greater 
than the collision time and the sudden approximation becomes accurate. 

Then we consider the sudden approximation in which u = and Eq. (J2J) 
can be calculated with the nuclear and Coulomb potentials to all orders 
giving 

g s ^ d (k, d; Im) = J d 2 r ± e~ i]l ^ (eT***^*"^ _ ^ ^ ( r±> k j ? (15) 
where 

/oo 
dt V eS (r, R(t)) = Qx (16) 

and V e ff is given by Eq. (|23|) . However as we mentioned above, in a first 
step we study only the effects of the Coulomb potential. Thus we call g^f dd 
the amplitude obtained from Eq. ([T5|) . by setting the nuclear potential equal 
to zero. Then the Coulomb amplitude in the sudden approximation to all 
orders can be written as 

^ udd (k,d;/m) = Jd 2 T ± e- ik ^(e- lQx -l)t lm (r ± ,k z ) 

= (film (k x + Qt ky, k z ) (fii m (k x , ky, k z ) (17) 

In the limit of very small Q Eq. (|17|) gives 

g s c uddl (k, d; Im) « Q-^-tim k y , k z ) (18) 

which is the sudden approximation restricted to first order and it agrees with 
the perturbation formula in the sudden limit, because mK\ (zu) = 1 and 
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wKq (zcr) = when w = 0. If the initial state wave function is approximated 
by its asymptotic form which is an Hankel function 



<piU^) = -i l Ci7ohi 1) (i'ror)Y lm (e,(f>), 7o r » 1, (19) 

where Cj is the asymptotic normalization constant and 70 = y/—2meo/h, 
then the general form of the initial state momentum distribution is given by 
the Fourier transform of Eq. (|T9~|) : 



(fc/7o)' 

(k 2 + 7o 2 : 



> ImW = AirCi ^ 1 ; u \ Y lm (k) (20) 



For the 2si/2 halo state of 11 Be it reads 

oo (fc) = 2^C 1 = 2yftC± (21) 
(k- + 75) K 2 

Where we have put K 2 = k 2 + 7 2 . By defining the dimensionless strength 
parameter \ — Q/K and using Eq. ()21jl. the sudden to all orders amplitude 
calculated explicitly up to second order in x reads 



^ dd (k,d;00) = 



2k x / Ak 2 x \ 2 



(22) 



In Appendix A, Eqs. ()29|) and (|30|) . we give the expansion of the amplitude 
up to second order in the full time dependent approach. Eq. (j22j) . can also be 
obtained from those equations in the limit zu — > 0. The strength parameter 
X represent the ratio of the classical Coulomb momentum transfer Q over the 
momentum K which is an average of the neutron final and initial momenta. 
If x is small, then a first order theory is accurate. From Eq. (J22j) one sees 
that if k x = 0, which happens for example for scattering at zero degrees, 
the amplitude gets contribution starting from the second order term. Thus 
the higher order terms are important for a proper description of the forward 
angle neutrons. 

The transition from the perturbative to the non-perturbative but all or- 
ders regime can be studied by plotting the strength parameter % vs the adi- 
abaticity parameter vj using the impact parameter as a variable. We show 
in Fig. Q the results obtained for various combinations of incident beam 
energy and neutron separation energy and different targets as indicated in 
the figure. The case of low incident energy clearly needs an exact treatment 
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Table 1: For several reactions discussed in this paper, we give the associated 
strong absorption radius R s and the critical impact parameter d crit . The 
maximum value of the impact parameter is always d = 120 fm. 



Target 


E inc (A.MeV) 


R s (fm) 


d crit (fm) 


u Be Be 


41 


6.0 


4.4 


Ti 


41 


8.2 


10.2 


Au 


41 


11.3 


19.4 


Pb 


72 


11.4 


19.8 


19 C Pb 


67 


12.0 


18.4 



of Coulomb breakup because at all impact parameters both the strength pa- 
rameter x as we U as the adiabaticity parameter w get values close to one. 
In the other cases instead, there is always one of the two limits which works 
well. For small separation energies (O.lMeV) w is very small and one can 
use the sudden approximation to all orders. On the other hand for large 
separation energies (5MeV) x is always small and the first order perturba- 
tion theory is accurate enough. In the cases we are discussing in this paper 
there is a smooth transition from one regime to the other and the transition 
occurs for impact parameters such that x ~ w which is satisfied in our case 

for d = d cr it = J Z err -/\/27o|eo|- This discussion and our formulation of the 
sudden approximation are very close in spirit to the work of Typel and Baur 
PU and to ref.[IE]. 

5 Results 

We discuss now the behavior of ■ An interesting observable which has 

been measured in a few experiments |2*H1 is the exclusive neutron angular 
distribution. The data of were taken at 41 A.MeV therefore it would be 
interesting to check if our model works reasonably well at such low energy. We 
show in Fig. (J2J) (a) and (b) angular distributions for the Coulomb breakup 
alone calculated for the reaction 11 Be+ 197 Au at 41 A.MeV at two impact 
parameters d = R s , and d — 50 fm. The solid line is the result obtained 
from gQ BM1 which is compared to the result from #cn M2 without nuclear 
interaction (big dots), the dotted line is the first order perturbation theory 
calculation, dashed line is the sudden to all orders while dot-dashed line is 
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the sudden to first order result. We see that at small impact parameters 
the sudden to first order and the sudden to all orders give coincident results, 
starting from an angle of about 10 deg. At d = 50 fm already the sudden to 
all orders and the sudden to first order give the same results at all angles. 
Thus it appears that higher order effects are important only at small neutron 
angles and at small impact parameters. Also in these situations the gQ BMl 
amplitude Eq. (fTTj) reduces to the first order perturbation theory amplitude 
Eq. (JT3j) . Eq. (|11|) can be considered correct only if the treatment of higher 
order terms by the sudden approximation, is accurate. This is indeed shown 
in Fig. © by the results obtained with the effective formula corrected to 
second order #cn M2 (big dots) calculated without nuclear potential. 

At angles smaller than 10 deg and small impact parameters the all orders 
calculation is different from the first order calculation. This is due to higher 
orders terms in the Coulomb field. Fig. © clarify this point: we represent 
the contribution of the first and the second order calculation normalized by 
their sum. It is clear that second order is very important at small neutron 
angles for small impact parameters whereas at large impact parameters, its 
effects are negligible. 

The angular distributions integrated over impact parameters are shown 
in Fig. (HJ) (a) and (b). The result with the g BBM1 amplitude in (a) is shown 
by the full line while the result with the Coulomb first order amplitude is 
the dotted line. In the case of pure Coulomb breakup, the first order ampli- 
tude is accurate enough down to 10 degrees, where higher order terms in the 
Coulomb field flatten the angular distribution. Such behavior seems indeed 
to be present in the experimental data which are shown on the right hand 
side figure together with the nuclear contribution. Our results explain why 
first order perturbation theory has been so successful in earlier studies of 
Coulomb breakup and they provide also a further justification of our pre- 
vious approach [llj. Note that for lighter targets the data do show a peak 
slightly shifted from zero degrees, thus reflecting the less important effect of 
higher order contributions [23]. In Fig. (J3J) (b) we give instead by the solid 
line the results from #cn M1 including Coulomb and nuclear potential. For 
completeness we give also by the dashed line the nuclear contribution alone 
and the experimental data from 

We present now results relative to energy distributions obtained for the 
reaction 11 Be( 208 Pb, 208 Pb) 10 Be+n at E- mc =72 A.MeV [23 . In Fig. © we give 
the results of the calculations obtained by the amplitude g™ 1 ! including 
both nuclear and Coulomb potential, Eq. (|12|) . by the solid line. The sudden 
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all-order amplitude (dot-dashed line) Eq. ()15j). the first order perturbation 
theory (dotted curve) Eq. (jHj) and the sudden first order (short dashed) 
Eq. ()18p for an impact parameter corresponding to the strong absorption 
radius d = R s and to d=20, 30 and 50 fm. The long dashed line gives 
only the nuclear contribution. At small impact parameters the two first 
order calculations, sudden and time dependent perturbation theory are very 
close, thus showing that the sudden approximation is valid at small impact 
parameters and therefore suggesting that Eq. (fT7|) would be accurate at small 
impact parameters to calculate the all orders amplitude. On the other hand 
we see that starting from d = 30 fm a new regime applies in which the two 
sudden calculations, to all orders or to first order give the same results. Then 
we conclude again that at high impact parameters the higher order effects 
can be neglected and first order perturbation theory applies. 

The main results of our new calculations are shown in Fig. (jUJ) (a) and 
(b) by the solid thick line. The effective formula gc^ M1 nas been integrated 
in the impact parameter range smaller and larger than d cr u = 19.8 fm re- 
spectively. These calculations indicate that for R s < d < 19.8 fm the results 
obtained with gQ^ M1 , solid line, are smaller than those obtained with first 
order perturbation theory g c ertl , thus showing higher order terms need to be 
considered. Since from the angular distributions shown in the previous fig- 
ures one can see that the effect of the higher order terms in Coulomb breakup 
alone is rather small, we suggest that the strong depletion shown by the peak 
of the energy distributions including Coulomb and nuclear breakup, comes 
mainly from the destructive interference effect already discussed in [11] . On 
the other hand for d > 19.8 fm we find that higher order effects are negligi- 
ble since using the effective formula or perturbation theory gives very little 
difference. We have checked that for 

d > 30 fm perturbation theory agrees exactly with the effective formula. 
Then we can conclude that at large d perturbation theory is valid. It is 
then reasonable to think that experimental data for neutron breakup could 
be analyzed by first order perturbation theory provided one could extract 
the contribution from impact parameters somewhat larger than d cr n. As we 
mentioned before, it is important to notice that the amplitude defined as 
9c , valid at all core-target impact parameters does not give rise to any 
divergences in the final integral over impact parameters. This is because 
the first order sudden term <7c uddl which contains the divergence is removed 
and substituted by the first order time dependent perturbation theory term 
gQ Tt 1 which does not diverge. The dashed line gives the nuclear contribution 
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alone and the dot-dashed line is the result of our previous method [TTj, re- 
calculated using relativistic kinematics for consistency with all other results 
of this paper. 

Finally the results of the full impact parameter integration are shown in 
Fig. (J7J) which gives the neutron final energy spectrum with respect to the core 
for breakup of n Be and 19 C on 208 Pb at 72 A.MeV and 67 A.MeV respectively. 
Experimental data are from [2U I2H1- Notations are as in Fig. (6). In the 
case of n Be the theoretical calculations have been multiplied by the known 
spectroscopic factor C 2 S = 0.77, while for 19 C we have used C 2 S = 0.65 and 
a neutron separation energy for the 2s state of 0.5 MeV. As expected, and 
already shown by other authors the effects of higher order terms are to reduce 
the peak cross section PQ. Analysis of the type presented in this section have 
been used to extract spectroscopic factors. 

6 Conclusions 

In this work we have presented an approximation scheme which allows cal- 
culating Coulomb and nuclear breakup and their coupling to all orders. It 
appears that higher order terms can give some effect only for heavy targets, 
at low core-target impact parameters and small neutron angles, and their 
effect is still noticeable after integrating over impact parameters. We have 
shown that higher order effects are accurately calculated by the sudden ap- 
proximation. The neutron angular distribution observable is well reproduced 
by the time dependent perturbation theory of Coulomb breakup. Once again 
we would like to stress the usefulness of measuring neutron angular distri- 
butions following breakup, as one of the best observables to clarify reaction 
mechanisms and to test theoretical models. 

The neutron-core relative energy spectrum after breakup shows a deple- 
tion of the peak value when higher order effects are included. It appears 
that this is mainly due to the destructive interference between nuclear and 
Coulomb breakup. Our conclusions are in agreement with recent numeri- 
cal solutions of the breakup problem Pfl 121 IE] and with new experimental 
data |2II] ■ In particular as higher order effects in the Coulomb potential are 
important only for small impact parameters (d < d C rit)> i n order to extract 
spectroscopic factors without any ambiguities, it appears very well suited 
to compare Coulomb first order perturbation theory with data containing 
contributions from d > d crit only j2Hj- 



12 



One neutron halo breakup is a simple reaction mechanism for which our 
approximation scheme allows to treat higher order terms in the Coulomb 
and the nuclear fields at the same time, including their couplings. This 
approximation scheme leads to simple expressions which can be generalized 
to proton breakup, where one needs to include the particle-target Coulomb 
potential and its quadrupole component. 

Aknowledgments 

We wish to thank Stefan Typel for discussions and comments and for 
communicating his results previuos to publication. 



A Eikonal perturbation theory for Coulomb 
breakup 

We give here some explicit formulae for the eikonal perturbation theory in 
the case of an electromagnetic excitation. The Coulomb field from a target 
nucleus can act both on the core and on the halo nucleus. Here we are 
only interested in the part that acts on their relative position and cause the 
breakup. For this reason, we subtract the part that acts on the position of 
the center of mass and we obtain [19j : 

V* = Z T e 2 ( , Zh , + , Zc , - (23) 
V|R + /3 2 r| |R-/3ir| |R| J K } 

where charges and masses are: core (A C ,Z C ), halo (A h ,Z h ), target (A T ,Z T ). 
We used also two ratios : /3i = Ah/Ap and = A c /Ap, with Ap = Ac + Ah. 
We develop this interaction in series and we take the case of a zero charge 
halo thus getting the dipole field 

V eS (r,t) = f(t)x + g(t)z (24) 

where 

fU\ ^csd Z eS vt 



(d 2 + (vtyf 2 ' (d? + ivtyf 2 ' 

and Z eS = p 1 Z P Z T e 2 . According to the eikonal formalism, the phase shift is 
expressed in the simple form : 

Xes{r,t) = J t V ee {r,t')dt' = F{t)x + G{t)z (25) 
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with 



F(t) = if /(*')*' = | (i 



Gffl = 1/^(0*' = !^ (26) 



where s = vt/d and Q = 2Z e g/ '{hvd). 

Introducing Eq. (|23)) into Eq. (jlj) leads to 



^ ert (k,d;/m) = -Jdte^j t 4> lm (k t ) (27) 

where k t = (A;.,, + F(t), k y , k z + G(t)). Using Eq. (J2T]) the breakup amplitude 
from an initial s-state becomes 

S r ( k,d;00) =4C,^ fdte- Mfe + F Wj + gWfe + gW] f . (28) 

J ([t + F(«)] 2 + fc, 2 + [*, + G(«)] 2 + 7 „ 2 ) 

We obtain the eikonal perturbation theory from an expansion in powers 
of the phase shift. Due to the simple form of the phase shift Xes = Qx, we ex- 
pand Eq. (|2*8j) with respect to the small quantities F{t)k x /K 2 and G(t)k z /K 2 
which are proportional to the parameter x = QjK used in Eq. (J2~2"j) . and ne- 
glect terms of order F(t) 2 and G(t) 2 . Thus the first order term reads : 

^ ertl (k,d;00) = -V^C i Q^(2k x K 1 (zj) + 2ik z K (zoj) (29) 

where one sees that the adiabaticity parameter w appears naturally. The 
second order amplitude is 

^ ert2 (k,d;00) = -^C^^^-Ak^K^-AiKKK^w) 

-m(k 2 + 2k x k z - k 2 )e-™) (30) 

The sudden approximation of these amplitudes are simply deduced as the 
limit w — > 0. 

In order to make a link with previous work 1J by other authors, we 
calculate explicitly the probability momentum distribution in the sudden 
approximation as: 

« = i/^W"Mf (3D 
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Then, we find for the first order and second order terms in the dipole field 
dP i( d ) _ 1 [ r]PJk 2 \n snddl (k till 2 - X (12) 

dP 22(d) _ 1 / i ^ 2b s U dd2 | 2_ V x 2 (15 - 10x 2 + 23x 4 ) (33) 



dk 8tt 3 J ^ 1 157T7o (1 + x 2 ) 6 

^13 (d) _ 1 f ^y^2gsuddlgsudd3 _ 32y 4 x 4 (5 — X 2 ) 



dk 8n 3 J " * * 15tt7o (1 + x 2 ) 

dPjjd) dPj 2 (d) dP^jd) Ay 4 x 2 (15 - 90x 2 + 39x 4 ) 
dk dk dk 15tt7o (1 + x 2 ) 6 1 ' 

with x = k/jo, y = Q^'/jo an d we used the asymptotic normalization 
constant C 2 = 2j of authors 0. 

Our first order result Eq. (|32|) is identical with the equivalent term Eq. (2.6) 
of ref . pQ . The second order dipole-dipole term of this work Eq. (|33|) is dif- 
ferent from Eq. (2.8) of 1 . The difference is due to the fact that in PQ the 
exact neutron-core scattering wave function 

i ikr 

/ (k)=e ik ' r - — . (36) 

7o — ik r 

was used, while in this paper the plane wave approximation has been applied, 
consistently with our hypothesis that the neutron-core final state interaction 
can be neglected (cf. Sec. 2 of [H]). We have checked that using Eq. (|3l)|) 
we would get the same result as in (I] but also that the difference with 
Eq. is negligible. To show this point we give in Fig. (JHJ) the sudden first 
order Coulomb breakup probability Eq. (j32)l after momentum integration as 
a function of the core-target impact parameter d (solid line) and the second 
order dipole-dipole term Eq. (J3l?j) calculated according to Eq. © (dashed 
line) and with the final plane wave function substituted by Eq. (j3Tlj) (short 
dashed line). It is clear that second order terms are rather small compared 
to the first order term, but also that the use of final plane waves is justified. 



B Fourier Transforms for the Coulomb po- 
tential 



J dte iuJt f(t) = QwK^w) (37) 
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dte iuJt g(t) = iQwK (zu) 



(38) 




Q to *r ( \ - n - 
— — I 2wKi{w) — t—zue 



4 V 2 



(39) 




) 



(40) 



(41) 



(42) 
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Figure 1: Perturbation parameter x versus adiabaticity parameter zu for 
impact parameters between R s and 110 fm and for the reaction 11 Be+Target. 
We vary the target (Pb, Be), the energy of the beam (72, 20 A.MeV) and 
the binding energy of the neutron (-0.1, -0.5, -5 MeV). Four regions where 
different theories should be valid are indicated: (a) sudden first order, (b) 
first order time dependent, (c) sudden to all order, (d) time dependent to all 
orders. 
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6 (deg) (deg) 

Figure 2: Neutron angular distributions for the reaction 11 Be+ 197 Au at 41 
A.MeV at d — R s and 50 fm. Solid line is the result of gQ BMl compared to 
#cn M2 calculated without nuclear potential (big dots). The dotted line is the 
first order perturbation theory calculation, dashed line is the sudden to all 
orders while dot-dashed line is the sudden to first order result. 
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Figure 3: Same reaction as Fig. (2), the lines are ratios of differential cross 
sections as a function of neutron angle. The solid and dotted lines are the 
ratios of the first and second order time dependent perturbation theory cross 
sections respectively divided by their sum. 
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Figure 4: Same reaction as Fig. (2) with the full d-integration. In both figures 
the dotted line is first order perturbation theory. On the left hand side figure, 
the solid line is the calculation with the effective amplitude including only 
the Coulomb potential, while on the right hand side the solid line includes 
both nuclear and Coulomb effects and the dashed line is the nuclear breakup 
alone. Data are from [2*3] . 



21 




Figure 5: Neutron-target energy distributions from n Be+ 208 Pb at 72 A.MeV 
for several impact parameters. Solid line is the result of gcN M \ the dotted 
line is the first order perturbation theory calculation, dashed line is the first 
order sudden calculation, the long dashed is the nuclear sudden to all orders 
and dot-dashed line is the Coulomb plus nuclear sudden to all orders result. 
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Figure 6: Same as Fig. (5) with integration over impact parameters smaller 
and larger than the critical value (cf. Table 1). The solid line is the #cn M1 
result, dotted is Coulomb first order perturbation theory, dashed line is nu- 
clear sudden to all orders and dotted-dash line is the previous result from 
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Figure 7: Comparison between calculations and experimental results for the 
reactions: (a) n Be+ 208 Pb at 72 A.MeV and (b) 19 C+ 208 Pb at 67 A.MeV 
(right). The solid line is the g^ M1 result, dotted is Coulomb first order 
perturbation theory, dashed line is nuclear sudden to all orders and dot- 
dashed line is previous results from [TT]. The spectroscopic factor used are 
C 2 S = 0.77 for n Be and C 2 S = 0.65 for 19 C. Data are from [2H I2S1 - 
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Figure 8: Comparison between the sudden first order Coulomb breakup prob- 
ability Eq. ()32|) after momentum integration, as a function of the core-target 
impact parameter (dashed line) and the second order dipole-dipole term cal- 
culated according to Eqs. (J2J) and (solid line) and with the final plane 
wave function substituted by Eq. (pffij) (dotted line). 
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